Low-lying excitations around a single vortex 
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A full quantum-mechanical treatment of the Bogoliubov-de Gennes equation for a single 
vortex in a d-wave superconductor is presented. First, we find low-energy states extended 
in four diagonal directions, which have no counterpart in a vortex of s-wave superconduc- 
tors. The four- fold symmetry is due to 'quantum effect', which is enhanced when is 
small. Second, for pp£, ~ 1, a peak with a large energy gap Eq ~ A is found in the density 
of states, which is due to the formation of the lowest bound states. 
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After a few years of controversy, d-wave nature of the high-T c superconductors is now 
well established fl], §], although superconductivity in the electron-doped Nd2_ x Ce a: Cu04 
appears to be of s-wave 0. Therefore it is important to understand the nature of vortex 
states in a d-wave superconductor 0, f|. An earlier analysis of the vortex state based on 
the Gor'kov equation shows that a square lattice of vortices tilted by 7r/4 from the a- axis is 
the most stable except in the immediate vicinty of T = T c or in a weak magnetic field ||. 
Such a square lattice of vortices, though distorted, has been seen by a small angle neutron 
scattering || and a scanning tunneling microscopy (STM) ]7| in YBa2Cu307_<5 (YBCO) 
at low temperatures (T < 10K). We believe that this distortion of the vortex lattice is due 
to the orthorhomibicity of the YBCO, although there are alternative interpretations based 
on the (d+s) admixture || ||. One of the most remarkable results in the STM experiment 
is that the vortex appears to have a circular symmetry as in an s-wave superconductor. 
It is in sharp contrast to earlier results obtained within the Eilenberger theory (a semi- 
classical theory of a superconductor) ||10|| , where a clear four-fold symmetry was obtained 
in the local density of states [11], [12|, [13|]. Further, at the center of the vortex, a peak 
with a large energy gap E ~ A was found in the local density of states, where A is 
the superconducting order parameter. Then the most natural interpretation is that this 
corresponds to the lowest bound state for a vortex in a d-wave superconductor analogous 
to the one predicted by Caroli, de Gennes and Matricon ||T4 |. 

In the previous study [15], in order to understand the results from the STM experiment, 
we have solved the Bogoliubov-de Gennes (B-dG) equation for a d-wave superconductor 
and obtained quasi-particle spectra around a single vortex. In the temperature region 
where the Ginzburg-Landau (GL) theory is valid, we found that the local density of states 
exhibits a circular symmetry and a peak with a large energy gap Eq ~ A is found in 
the local density of states at the center of the vortex, which is consistent with the STM 



experiment. In ref. |J5], it is crucial to set pp£, ~ 1 for YBCO, where pp and £ are the 
Fermi momentum and the coherence length respectively. The value of pp^ is obtained 
by an approximate formula for the lowest bound state E = A/(7ipp£) |12], |IB| . This is 
also consistent with the chemical potential of YBCO deduced from the analysis of the 
spin gap seen in an inelastic neutron scattering experiment from monocrystals of YBCO 
|17], In the above analysis, however, we neglected the noncommutability between k and 
x ('quantum effect where k and x are the quasi-particle momentum and coordinate 



respectively, and the local density of states has a perfect circular symmetry except when the 
mixing of an s-wave component occurs [pOR . The correction is 0(l/pp^) and irrelevant at 
least in the study of systems with a long cohernce length e.g. the superconducting phases 
of the heavy-ferimon systems ~ 10) and the 3 He superfluidity ~ 100), but may 
have a serious influence in the study of the high-T c superconductors, where pp£, ~ 1 as is 
discussed above. 

In this paper, a full quantum-mechanical treatment of the B-dG equation for a d-wave 
superconductor is reported, where the 'quantum effect' is taken into account. As shown 
below, the four-fold symmetry appears in the local density of states. Similar four-fold 
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symmetry was obtained in the previous studies ]TT], [L2|, [Tj|. But it should be noted that 
the four-fold symmetry discussed here has totally different origin from that obtained in 
the earlier studies. The B-dG equation for a d-wave superconductor is given by 



{- — (V-ieA(K)) 2 -ii}u n (K) 
-{d x {A(x)d x ) - d y (A(x)d y )}v n (x.) = e n w„(x), 
-{-^(V + zeA(x)) 2 -/i}^(x) 
-{d x (A(x)*d x ) - d y (A(x)*d y )}u n (x) = e n v n (x.), 



where u n (x) and u n (x) are the quasi-particle amplitudes, A(x) is the pair potential, A(x) 
is the vector potential which is neglected assuming H<^iH C 2, and fi is the chemical potential 
which is identified with the Fermi energy The parameters are set as 2m£ 2 A = 2.82 and 
R/£ = 30, where £ is the coherence length and R is radius of the system and the pair 
potential A(x) is given by 

A(x) = Atanh(r/£)e^, (0.1) 

where r and are defined by x = (r cos(0), r sin(0)). This form of pair potential is obtained 
by the Ginzburg-Landau (GL) theory for a d-wave superconductor [T|, ^0f and applicable 
at not too low temparatures (an estimate of the temperature region where the GL theory 
is valid is [0.5T C ,T C ] for an s-wave superconductor, when one consider the quasi-particle 
spectra around a single vortex). 

In order to solve the B-dG equation numerically, it is convenient to expand the quasi- 
particle amplitudes w n (x) and t> n (x) as 



oo oo 

u n (x) = XX^i>l*l( r ) ex PW)> 

l=— ooj=l 
oo oo 

v n(k) = &™>'jA|z-i|( r ) exp(i(Z - !)</>)• 



l=—ooj=l 



Here ipiJx) = 7^=^-r , — ^Juio-ivxl R) (J v (x) is the Bessel function), ola v is the 7-th 



positive zero point of J v (x) and R is the radius of the system. In the previous study |L5 
in which the 'quantum effect' is neglected, the B-dG equation decouples for m„,zj's and 
v n,i,j's with different I. On the other hand, when the 'quantum effect' is taken into account, 
Wn^j's and v n j/s with different I couple, which plays an important role when p F £ is small. 
Because of the coupling, the number of basis we can use for numerical diagonalization is 
small compared to the previous study |L5|]. However, it is sufficient for the understanding 



of the qualitative aspects. 

At first, consider the density of states J2i${E — Ei) as a function of E/A, where 
Pf£, = 1.33. As is seen in Fig. 1, there is a peak with a large energy gap Eq ~ A. This 
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corresponds to the lowest bound state. The peak has a width due to the internal degree of 
freedom in the k space. These are consistent with the previous results |L5| qualitatively, 
where the 'quantum effect' is neglected. 

Next consider a local density of states in a superconductor, which is the quantity of 
interest for comparison to STM experiments and given by, 

N(E,x) =£[ K(x)| 2 S(E-e n ) 

n 

+ K(x)| 2 5(E + e n )}. 

In Fig. 2 and 3, / a35A dE N(E, r, 4>) is plotted and they show a clear four-fold symmetry 
in the local density of states. When p F £ is changed from 1.33 to 2.00, the four-fold 
symmetry is suppressed (see Fig. 3) and the local density of states becomes circular. 
This supports the idea that the four- fold symmetry is due to the 'quantum effect'. We 
stress that these low-energy states extended in four diagonal directions are paticular to 
d-wave superconductivity. Therefore 1) these states should give rise to the zero-energy 



density of states proportional to \JB as discussed by Volovik and others ^TJ f4j, 2) they 
are the most likely the origin of the large flux flow resitivity in YBCO observed recently 
by Doettinger et al. |2^] and, 3) when a square lattice of vortices tilted by 7r/4 from 



the a-axis is formed, the quasi-paricle can move from one vortex to the other through 
these low-energy states extended in four diagonal directions, which should give rise to a 
cohesive energy guaranteeing the stability of the square lattice. The clarification of the 
quasi-particle spectrum in a vortex lattice and, in paticular, the tilted square lattice is of 
immediate interest. 

In conclusion, we have investigated the Bogoliubov-de Gennes equation for a d-wave 
superconductor, where the noncommutability between k and x ('quantum effect') is taken 
into account. We found a peak with a large energy gap E ~ A in the density of states, 
which is consistent with the previous results. We found low-energy states extended in four 
diagonal directions, which is due to the 'quantum effect'. The low-energy states have no 
counterpart in a vortex of s-wave superconductors. It is natural to consider that these 
low-energy states cause directional attractive forces between vortices. It is possible that, 
due to the directional attractive force, a square lattice of vortices becomes stable in some 
parameter region. Another scenario for a square lattice of vortices is proposed in ref. ||, 



where the higher-order correction in the Ginzburg-Landau theory |23| plays an essential 
role. In this paper, we do not consider the effect of the higher-order correction. The 
higher-order correction causes the four-fold symmetry in the pair potential. We consider 
that, in low temperatures, it is needed to take into account both the 'quantum effect' and 
the higher-order correction in the GL theory self-consistently, and more detailed study is 
left as a future problem. 

This work is in part supported by National Science Foundation under grant number 
DMR95-31720, ISI Foundation at Torino, Italy and timely short term research fellowship 
of Japan Society of Promotion of Science. 
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Figure Captions 

Fig. 1: Ei 6{E - Ei) as a function of E/A, where k F £ = 1.33. 

Fig. 2: Jo' 35A dE N(E,r,<f>), where k F £ = 1.33 for a) r/f = 3.0, b) r/f = 9.0, c) 
r/f = 15.0 and d) r/f = 27.0. 

Fig. 3: Jo- 35A dE N(E,r,4>), where a) £; F £ = 1.33 and r/f = 9.0 and b) k F £ = 2.0 and 
r/i = 9.0. 
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Fig. 2 



— i 1 1 1 1 1 r 

o o o 

o oo o© oo © 

o © © © © © © ©_ 

© ©© ©© ©© © 

© ©© ©© ©© © 

ooooooo© 

oo o© ©© ©© 

©o 000 ©© ©0°°°©© ©0 000 ©© ©o 000 ©* 

I I I I I I L 

1 2 3 4 5 6,7 



Fig. 3 a) 



N 




Fig. 3 b) 



